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INTRODUCTION 


Researchers have provided forest managers and their consultants with several 
computer programs for analyzing investments in forestry (Hall 1962, Row 1963, 
Schweitzer et al. 1967, and Forster 1968). These computer programs have been 
modified in many cases to apply to situations not originally envisioned by the 
program developers (e.g., see Marty et al. 1966, Green and Alley 1967, Wikstrom 
and Alley 1968). These programs are particularly useful if the situation being 
analyzed fits one of the options available in the program and if the manager's 
investment criterion is handled by the program. However, most of the available 
computer programs are cpomoned for specific problems and, therefore, are rela- 
tively complex to use.!/ In addition, only two investment criteria are available 
in the programs; namely, internal rate of return and ,present net worth. 


This paper is meant to fulfill the following objectives. First, to announce 
the availability of the IVST computer program to prospective users. Second, to 
discuss the applicability of this program to managerial problems. Finally, to 
outline the operating characteristics of the program so analysts can use it. 


PROBLEMS WHICH MAY BE ANALYZED WITH IVST 


IVST, in common with all other investment analysis computer programs, is 
designed to help the decisionmaker decide whether or not to pursue certain pro- 
ductive activities or alternatives. There are many types of investment decisions, 
including: (1) whether a certain kind of productive activity should be pursued 
at all, (2) the intensity at which a productive enterprise should be operated, 

(3) timing of the productive activity in relation to other activities, (4) 
sequence of the productive activity in relation to other activities, (5) the 
location at which the productive activity takes place, and (6) the best way to 
accomplish the productive activity--i.e., technology. Depending upon the way 

the alternatives are formulated, any decision may be included in an IVST analysis. 


The type of decision to be analyzed dictates the kinds of data which are 
required. Here it is presumed that data will be expressed in economic units, 
e.g., dollars. Data will consist, therefore, of costs and returns tagged with 
a time dimension--date of occurrence. Realistically, the costs and returns must 
be regarded as managerial expectations, since an investment analysis of histori- 
cal occurrences provides few guidelines to the future. 


However, data alone are not sufficient to guide managers in their decision- 
making function. Also required is a decision rule or criterion, which simply 
stated is the manager's systematic way of evaluating how well alternative strat- 
egies provide for fulfilling the firm's goals. 


Now, whereas the type of investment is completely open when IVST is used, 
the choice of criterion is not. Three alternative investment criteria are 
available in IVST: internal rate of return, present net worth, and benefit-cost 


| the rate-of-return computer programs devised by Hall (1962) and Forster 
(1968) are exceptions, since they are general and presume no particular problem 
context. 


ratio.2/ It is the job of the manager to choose the criterion which will best 
provide for fulfillment of the firm's goals. If one of the three criteria avail- 
able in IVST will fill this bill, then this computer program will be useful to 
the manager. On the other hand, if the manager selects an investment criterion 
that is not available in IVST, then the program will not be useful and, indeed, 
its use would be harmful, since misallocated resources would result. 


IVST is a generalized computer program in that no particular problem context 
is presumed by its structure and output formats. The program is simple to use 
in that input requirements are minimal for the types of problems which may be 
analyzed. It is assumed that all preliminary data preparation is handled exter- 
nally. Input data consist of a time-tagged stream of costs and returns and 
certain identifying information. A data preparation computer program will be 
required if data are in the form of quantities and prices. Large problems may 
be analyzed since data storage requirements are low per alternative analyzed. 
Finally, as pointed out above, IVST permits the analyst to select for examination 
one or more of three optional investment criteria: present net worth, internal 
rate of return, and the benefit-cost ratio. In addition, the structure of IVST 
is such as to permit other criteria subroutines to be added if desired by the 
analyst. 


IVST was derived from two computer programs developed to analyze data for 
the Douglas-fir Supply Study, a cooperative study by the Pacific Northwest Region, 
Forest Service, and this Station. One of the computer programs was based on the 
internal rate-of-return criterion, the other on the present net worth and benefit- 
cost ratio criteria. These computer programs were modified to become options 
within IVST. 


STRUCTURE OF IVST 


IVST consists of a main program and two subroutines. The main program reads 
all data and provides for the printing of the first report, which consists of 
the input data. Subroutine BCA calculates the benefit-cost ratio and present net 
worth for the range of interest rates selected by the user. Subroutine IROR 
calculates the internal rate of return.3/ All results of the analyses are printed 
under control of the subroutines. Either BCA and/or IROR may be selected for any 
problem by the analyst. 


Both BCA and IROR are closed or linked subroutines, i.e., a subroutine not 
stored in the path of the main routine, IVST. A closed subroutine is entered by 
a jump operation in the main routine and return to the main routine is provided 
at the end of the operation. Because of this configuration, it is possible to 


2/ The reader interested in the theoretical aspects of these economic cri- 
teria should consult Eckstein (1958, pp. 70-79) on the benefit-cost ratio, and 
Hirshleifer (1958, pp. 329-352) on the internal rate of return and present net 
worth criteria. A comparison of management guides developed by application of 
these criteria to a forestry problem may be found in Webster (1965). 


3/ Development of IROR began with a major modification of Row's rate-of- 
return program (Row 1963). The author wishes to acknowledge the programing 
assistance of Edgel E. Skinner of the Station's statistical and data processing 
services staff in this initial effort. 
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make use of the subroutines independently of the main routine. That is, the 
prospective user may design his own main routine and simply call either BCA 

and/or IROR when needed. It may be desirable to do this when input data are not 

in proper form for the present main routine. Prospective users are urged to 
consult a computer programer if the coding of a new main routine appears necessary. 


If the user wishes output for all three criteria options, then a total of 
four reports will be printed. As noted above, the first report consists of the 
input data. In the second report, the benefit-cost ratio and present net worth 
for every interest rate within the selected range will be listed for each alter- 
native. The third report consists of certain messages generated within the IROR 
subroutine that will assist the analyst in interpreting results of the internal 
rate-of-return analysis. Finally, the fourth report lists present net worth for 
each interest rate within the selected range for each alternative until the in- 
ternal rate of return is found. The internal rate of return is listed as the 
last interest rate for each alternative. After the last report is printed, con- 
trol is switched back from the subroutines to the main routine and an end-of- 
problem card is read. At that time, the user can either terminate the program 
or solve another problem. 


The Main Program IVST 
Function 


The function of the main program is to read a stream of periodic annual 
costs, a stream of periodic annual returns, and associated data; prepare a report 
of the input data; generate annual costs and annual returns if data are for a 
period other than 1 year; and call in the criteria subroutines desired by the 
user. IVST is designed to accept periodic annual data for whatever period length 
the user desires. However, the same length period must apply to each cost and 
return entered for that particular problem. Of course, if period length differs 
from one alternative to another, each alternative may be run as a separate prob- 
lem. If period length varies within an alternative, however, it is mandatory 
that period length be set at 1 year and that annual data be read. 


The advantages of having the capability of reading periodic annual costs 
and returns, as opposed to annual inputs, include savings of time, card costs, 
keypunching costs, the minimization of input error, and convenience. 

IVST permits initial costs and returns to be read when a period length of 


1 year is designated. Initial costs and returns are those costs and returns that 
occur in year zero of tne investment series. 


IVST Inputs 


There are 17 specific types of inputs to the IVST program. Formats of in- 
put cards are shown in table 1 of the Appendix. The inputs are: 


1. Alphanumeric study identification. 
2. Beginning interest rate of the range to be examined, decimal. 


3. Final or ending interest rate of the range to be examined, decimal. 


Uo 


Wo 


Uso 
14. 
Us 
16. 
We 


Increment in interest rate, decimal. 

Number of alternatives in problem. 

Indicator of an initial cost or return (i.e., a cost or return occur- 
ring in year 0). If initial costs and returns are used, then period 
length must be 1 year (i.e., annual data must be read). 

Alphanumeric identification of alternatives. 


Problem number. 


Series type indicator (either terminable or perpetual). 


Investment criterion indicator (either internal rate of return or bene- 


fit-cost ratio and present net worth, or all three). 


Length of period over which periodic costs and returns apply, years. 
Annual data may be read by setting period length equal to 1. 


Maximum number of years in any alternative, i.e., length of longest 
investment series of any alternative in the problem. 


Number of years in each alternative. 

Alphanumeric problem identification. 

Periodic annual costs for each alternative, dollars. 
Periodic annual returns for each alternative, dollars. 


End of problem indicator. 


IVST Main Program Outputs 


Input items 1, 7, 8, 9, 11, 14, 15, and 16 constitute the first report of 
IVST (see fig. 1). Data for five alternatives are printed on each page. 


INVESTMENT ANALYSIS== RESEARCH PAPER SAMPLE PROBLEM 
PROBLEM NOD 1 TEST CATA 
TERMINABLE SERIES PERIGD= 5 YEARS 
MANAGEMENT SLTERNATIVE 
ONE TwC THREE FOUR FIVE 
PERIOD ANNUAL ANNUAL ANNUAL ANNUAL ANNUAL ANNUAL ANNUAL ANNUAL ANNUAL ANNUAL 
CccsT RETURN Cast RETURN CCST RETURN ccst RETURN CGST RETURN 
a 15 12 2266425 104254 32 34 33 Si 65 12 
2 23 65 2z92€CC 29237 45 52 - 22 $2 21 23 
3 65 42 1Z7S57C =T4550 es 12 44 $3 45 S2 
4 34 €3 ~21EEESS 1562246 7é 65 55 S4 23 45 
5 ie) a 2214230 1584476 Sa 8S 6€é S5 57 73 
é 0 ) -ZC21CS4 LE07ES4 14 35 88 sé 84 65 
7 0 ce) c Cc és 67 77 37 52 98 
8 0 i} c 9 23 85 22 se il €3 
Ss 0 i) Cc ie] 15 64 li 31 32 €5 
10 0 0 c c 24 13 22 a2 ll 32 
ll Q 9 0 0 Ss? 25 44 34 44 12 
12 QO ce} c c é2 35 66 25 55 45 
13 QO 0 c 0 54 53 bs 26 22 7 
14 ie) a} Ly} c 13 $2 hb) 37 33 $6 
15 0 0 Q ¢ 2 45 33 38 65 32 
16 re) ) c c c Q 55 3s 47 58 
1? c ie} ie} c c 0 45 22 0 a 
18 ie} 0 Q c C ° 23 45 2 
1s 0 3} c c c Q 15 €8 3 
20 0 ce} ¢ c Cc Q 72 2 2 


a 


Figure 1.--First report of the IVST computer program showing input data. 


Subroutine BCA 


Function 


The function of subroutine BCA is to calculate present net worth and benefit- 
cost ratio for each alternative for the range of interest rates selected by the 
analyst. 

Criterion 
Present net worth of a terminable series of costs and returns (dollars) is: 


R C 


Price emia en The ee ea (1) 
n t=O Qe t=O (ieee ye 
where 
s : 4/ 
Cz = cost in year t, dollars— 
R, = return or benefit in year t, dollars“! 
t = discount rate, decimal 
n = number of years in the investment series 


Present net worth of a perpetual series of costs. and returns is: 


= PVW. . (1 + i)” 2 
ne acs ea eee (2) 


In contrast, the benefit-cost ratio is: 


a (3) 


Equation 3 can be used for both terminable and perpetual series because the in- 
finite series multiplier (the last term in equation 2) cancels out of equation 3. 


Subroutine BCA does not assume that discounted returns (or benefits) or dis- 
counted costs are positive. It is likely that both negative returns and negative 
costs will occur if the analyst is working with marginal values (i.e., differ- 
ences from some base). 


4/ Any asset salvage values (revenues) or replacement costs should be in- 
cluded in the last term of the series. 


Four conditions are possible: 


Condition Return (or benefit) Cost 
1 + + 
2 - + 
3 + s 
4 a : 


Since the benefit-cost ratio must always be positive, it was necessary to 
include rules for conditions 2, 3, and 4 within the BCA subroutine. In conditions 
2 and 3, the B/C is set equal to zero. In condition 2, the reasoning is that a 
negative return is equivalent to a positive cost, but returns are equal to zero, 
hence B/C = 0/+ = 0. Likewise, in condition 3 a negative cost is equivalent to 
a positive return, hence cost is equal to zero. Therefore, B/C = +/0, which is 
undefined. These two conditions can be easily discerned since condition 2 will 
always be associated with a negative present net worth, whereas present net worth 
will always be positive in the case of condition 3. 


In condition 4, negative returns are equivalent to positive costs and nega- 
tive costs are equivalent to positive returns, hence the benefit-cost ratio is 
set equal to its inverse, the ratio of costs to benefits (C/B). 


In addition, it is possible that either discounted costs or discounted re- 
turns are equal to zero. In such a case, the BCA subroutine will set B/C equal 
to zero. 


BCA Outputs 


Subroutine BCA provides for printing the present net worth and benefit-cost 
ratio for each alternative for the entire range in interest rates, together with 
identification (see fig. 2) .2/ Results for five alternatives are printed on each 
page. After results are printed for all alternatives, control is returned to 
the main program. 


Subroutine IROR 


Function 


The function of subroutine IROR is to calculate the internal rate of return 
for each alternative within a range of starting interest rates selected by the 
user, using the annual costs and annual returns generated by the IVST main pro- 
gram. 


5/ The discounted costs and discounted returns may also be printed if the 
user so desires. The slight modification of the BCA subroutine which will pro- 
vide this output is described in the Appendix, page 17. 


PENEFIT=CCST ANALYSIS=~ RESEARCH PAPER SAMPLE PRCBLEM 
PRCBLEM NOe 1 TEST DATA 


BENEFIT-COST RATIO (@/C) AT ALTERNATIVE RATES CF INTEREST 
PRESENT NET WCRTH (PNW) AT ALTERNATIVE RATES OF INTEREST 


MANAGEMENT ALTERNATIVE 


ONE TRC THREE FOUR FIVE 
RATE >PNW B/C* *¥PNW B/C* *PNW B/C* *PNW B/C* *PNW B/C* 
De 0e325CE 03 1047 Co27G3E Ce Qo 924550E 03 1513 Oo1l155E 04 1,24 62935CE 03 1.28 
0259 00303ZE C3 1047 0.2202E 0& Oo 0e2284E 03 1.12 Ooll71E 04 1.31 Oe71L02E 03 1.26 
1.00 Co2832E 03 1.47 Col7ESE 08 12.19 022322E 03 1.10 Ooll59E 04 1.38 025356E 03 1323 
1,50 0.264EE C3 1.4€ O,13E5E ce 5ell 9.15e8aE 93 1,08 De1l39E 04 1245 Oo3991E 93 1.29 
2090 Oe247EE C3 104€ Co1CS4E Ce 3015 JslC26E 93 41096 OsLOSVE 04 1353 0,2915E 93 1,17 
2250 90232z2E 33 104€ Oo TE57E C7 2023 925956E 02 1,04 O.1045E 04 1.61 Oe2963E C3 1.14 
3200 Ceo217eE 93 1045 Co515CE 97 1.070 022650E 02 1,92 029984E 03 1.68 521383E 3 1.11 
3050 Qe204EE 03 1045 Oe2S71E 07 1.36 Jell@6E O1 1.99 O951CE 03 1.75 Oe8375E Uz 1.07 
4.00 0e1923E 63 1545 Co1l97e8E C7 1312 -Gel810E 22 0098 929048E 03 1,81 053972E 02 1,04 
4e5C UealBCSe 03 1044 “005625E O€é C394 —002268E 02 0.97 02.8695E 03 1.88 Oe4915E O1 1,0) 
5200 Qol7C4E C2 1244 “Ool1SE3ZE C7 Ooe@l —004355E 02 0395 0e8184E 03 1,93 —0.2506E U2 0.97 
505 DelENEE 03 1044 -9.3211E 07 0073 7925152E 92 0094 Qo7789E 03 1.95 —U2.4884E 02 0094 
6000 Del516E 03 1.044 -S,4270E Q7 Qebl ~Oe57i7E 02 6293 Oo7419E 03 2004 -0,6835E v2 9.91 
6259 901432E 03 1,43 —“JoS1G1lE 07 0.54 ~TeE19N1E 92 0092 0.7975E 03 2009 —008435E 02 0088 
7200 Qol354E 02 12423 ~CoS5SE2E 07 C048 ~9e633S5E 92 0091 006754E 03 2014 —-UeITSO1LE 02 0,85 
7250 Oel281E C3 1.42 ~3,E€62SE C7 0243 ~IJ26464E 02 0299 006457E 03 2018 ~CelI85E 03 0.82 
8.0C Cel213E 03 1.42 -O,71SEE 07 0035 -2264S5SE 02 0.89 Oo6180E U3 2022 ~Jell75E 03 0279 
8o5C QolL5CE 03 1.42 -0,7675E 07 0035 -006463E 32 0239 0,5923E 03 2026 -031249E 023 O.77 
9200 O.10SCE C2 1042 ~CoEQ7TeE C7? 0032 -026373E 02 Oe 88 0.5685E 03 2025 ~OQeol31CE 33 0274 
9250 Os.l025E 93 1242 -Co&412E 9 0030 -O2.6241E 92 0,88 005462E 93 2033 ~Gel360E 03 0272 
10,0C 009834E CZ 1.42 -25868SE C7 Co27 -OQ.6076E 0 0088 025255E 03 2236 -Os1401E 23 0.70 


Figure 2.--Second report of the IVST computer program showing results of 
use of the present net worth and benefit-cost ratio criteria. 


Criterion 

The internal rate of return (z) is defined as that interest rate which 
equates discounted returns to discounted costs, i.e., the rate that results in a 
present net worth of zero. If we use the same notation as above, then the inter- 


nal rate of return is defined as z when the following condition is satisfied: 


n Rt n Ct f 
ip pees 20 tA 
tzo (1+ 2)t a (+ ae for a terminable series (4) 


or 


n Ry (letaa) 5 ee Ges : (eee) (5) 


a6 (ee ae (Pe ae a1 Hees pene ya Rl = Be aa 


for a perpetual series. (1 +7)” 

It should be noted that infinite series multiplier |(1 +72)” - 1] cancels 
in equation 5, leaving it equivalent to equation 4. Therefore, the internal rate 
of return for a perpetual series is identical to the internal rate for a termi- 
nable series. However, discounted costs and returns will not be identical for 
both series. 


At this point it may be helpful to indicate how the three investment criteria 
relate to one another. When the discount rate is equal to the internal rate and 
hence discounted costs are equal to discounted returns (present net worth is equal 
to zero), then the benefit-cost ratio is one (compare equations 1, 3, and 4). 


Except for Forster's program (1968), all of the rate-of-return computer pro- 
grams cited earlier use unmodified iteration to isolate the internal rate of 
return. In contrast, in IVST present net worth is iterated at interest rate in- 
tervals selected by the analyst until a change of sign is found. At that point, 


the bisection method is applied until the internal rate of return is found with- 
in the interval at which the sign change occurs.6/ 7, 


The power of this technique can best be realized by selecting a wide range 
of interest rates and a large value for the rate increment. In this way, con- 
vergence to the internal rate of return occurs quickly as large bisection inter- 
vals are used. 


However, if one desires to calculate present net worth and the benefit-cost 
ratio as well as internal rate of return, the advantage of selecting a wide range 
of interest rates and a large value for the rate increment may be lost. The 
analyst will generally wish to use the BCA subroutine with a small value for the 
rate increment. This in turn may necessitate investigating a narrow range of 
interest rates (see program limitation 8, page 14). Because interest rates to 
be used in both subroutines are determined in the main program, they must repre- 
sent a compromise if the two subroutines are to be used together. Therefore, 
in cases such as this, the analyst may prefer to use the subroutines sequentially. 
This can be accomplished by setting up two problems for a single computer run, 
both with identical data except for interest rate range and increment. One of 
the problems would be run with the IROR subroutine with a wide range of interest 
rates and large increment, and the other problem would be run with the BCA sub- 
routine with a narrow range of interest rates and a small increment. 


An example which illustrates the bisection method utilized in IROR is shown 
in figure 3. In this example a range of interest rates of 0 to 10 percent at 
increments of 5 percent were specified. The program calculated the present net 
worth (PNW) at a rate of O percent and found it to be $455. Next PNW was calcu- 
lated at 5 percent and found to be -$43. At this point, the change in sign 
initiates the bisection procedure. The program calculates PNW halfway between 
0 percent and 5 percent, that is at 2.50 percent. Here PNW is positive ($60), 
so the program tries the rate midway between 2.50 percent and 5 percent, or 3.75 
percent. PNW at that rate is negative (-$9), so it is apparent that the internal 
rate lies between 2.50 percent and 3.75 percent. This procedure is continued 
until the internal rate of return is found at 3.52 percent on the sixth bisection 


(see fig. 3). 
IROR Outputs 


Three types of messages may be aes prior to the printing of tables 
showing the internal rate of return (see fig. 4). 


6/ For a discussion of the use of the bisection method for finding a real 
zero, see Hamming (1962, p. 352). The bisection method can result in both de- 
creased computation time and increased accuracy. As Hamming points out, "Since 
each step halves the interval in which the zero lies, ten steps will reduce the 
interval by a factor of about 1,000; 20 steps, 1,000,000; etc. Thus the method, 
which assumes only continuity and the ability to evaluate the function at any 
point, is fairly effective." 


Li} The recent rate-of-return computer program develcped by Forster (1968) 
also utilizes the bisection method, which he terms "interval halving." 
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Figure 3.--An example of how the IROR subroutine converges to the internal 
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DETERMINING RATE OF RETURN-— RESEARCH FAPER SAMPLE PROBLEM 
PRCBLEM NOs 1 TEST DATA 
MESSAGES 
INTERNAL RATE CF RETURN CCES NOT FALL WITHIN SELECTED RANGE FCR ALTERNATIVE 1 


PISECTION LIMIT EXCEECEC IN ALTERNATIVE ThC | 
CALCULATED FRESENT NET wWCRTH COES NCT FALL BETWEEN PLUS ANC FINUS 1200 COLLARS | 


PNW AS A FUNCTICN CF INTEREST RATE IS NCT TENDING TOWARD ZERO WITHIN A PCRTICN CF ITS RANGE FOR ALTERNATIVE 4 
Figure 4.--Third report of the IVST computer program indicating messages 
generated by the IROR subroutine. 


If the internal rate of return is not found within the range of interest 
rates selected by the analyst, the message 


INTERNAL RATE OF RETURN DOES NOT FALL WITHIN SELECTED RANGE FOR 
ALTERNATIVE X 


will be generated, with the alternative number filled in in place of X. 


The second message that may be generated informs the user of the precision 
levels that apply to the internal rates of return. Precision level is indicated 
only when bisection is used. If the internal rate is found on iteration, a zero | 
will appear in the PNW column and the precision level is to the nearest $0.01. | 


Three precision levels are built into the bisection phase of the IROR sub- 
routine, $1, $10, and $100 (i.e., rate-of-return equation identity may be correct 
either to the nearest $1, $10, or $100). Up to 100 bisections may occur within 
each of these precision levels. Whenever an internal rate fails to be determined 
at one of these precision levels, a message is generated and the next lower level 
of precision is tried. The message is of the form 


BISECTION LIMIT EXCEEDED IN ALTERNATIVE X. CALCULATED PRESENT 
NET WORTH DOES NOT FALL BETWEEN PLUS AND MINUS Y DOLLARS 


with the alternative number filled in by IROR in place of X, and the precision 
level in place of Y. 


A third message will be generated whenever present net worth as a function 
of interest rate does not converge toward zero. In this case, the form of the 
message is 


PNW AS A FUNCTION OF INTEREST RATE IS NOT TENDING TOWARD ZERO 
WITHIN A PORTION OF ITS RANGE FOR ALTERNATIVE X 


Occurrence of this message does not necessarily mean that the internal rate 
of return does not fall within the range of interest rates selected by the user. 
However, IROR will discontinue searching for the internal rate whenever the above 
condition is encountered. When this condition occurs, the analyst is advised to 
use the BCA subroutine to generate sufficient points to indicate the form of the 
present net worth function. If the above condition is caused by portions of the 
function which are not monotonic, it is likely that the interest rate range can 
be modified so as to exclude such portions, thereby permitting IROR to converge 
to the internal rate of return. 
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The behavior of IROR in various situations is illustrated in figure 5, 
wherein 10 representative types of present net worth functions are depicted .8/ 
The types shown in figure 5 by no means exhaust the possibilities. In each graph, 
present net worth (dollars) is plotted over the range of interest rates selected 
for examination. IROR will treat each type of function as follows: 


Type 1.--Present net worth is a monotonic decreasing function of interest 
rate. The internal rate of return falls within the range selected and will be 
determined by IROR. 


Type 2.--Present net worth is a monotonic increasing function of interest 
rate. The internal rate of return falls within the range selected and will be 
determined by IROR. 


Type &.--Present net worth is a monotonic decreasing function of interest 
rate. However, the internal rate of return is located outside the range selected. 


Type 4.--Present net worth is a monotonic increasing function of interest 
rate. However, the internal rate of return is located outside the range selected. 


In the case of types 3 and 4, IROR will calculate present net worth for the 
full range of interest rates and then generate the message 


INTERNAL RATE OF RETURN DOES NOT FALL WITHIN SELECTED RANGE FOR 
ALTERNATIVE X 


Type 5.--Present net worth is a monotonic decreasing function of interest 
rate. The internal rate of return is located outside the range selected, and 
present net worth is not converging toward zero as interest rate increases. 


Type 6.--Present net worth is a monotonic increasing function of interest 
rate. The internal rate of return is located outside the range selected, and 
present net worth is not converging toward zero as interest rate increases. 


In the case of types 5 and 6, IROR will calculate two iterated points along 
that segment and then print the message 


PNW AS A FUNCTION OF INTEREST RATE IS NOT TENDING TOWARD ZERO 
WITHIN A PORTION OF ITS RANGE FOR ALTERNATIVE X 


Types 7 and 8.--Present net worth is not a monotonic function of interest 
rate. Although the internal rate of return lies within the range of rates 
selected, IROR will not determine it. Rather, these situations will be handled 
by IROR exactly as types 5 and 6. If these types of functions are suspected, 
the analyst should examine PNW over a wide range of interest rates by use of BCA 
subroutine. After examining BCA outputs, it is possible to reset the range in 
interest rates to avoid the program stop, thereby permitting the IROR subroutine 
to converge to the internal rate of return. Also, it should be pointed out that 


8/ Samuelson (1937, p. 475) indicates that some investment opportunities 
may have no real internal rate of return (i.e., the present net worth equation 
has only imaginary roots). Also, he mentions the possibility of multiple internal 
rates, 1.e., the present net worth function may equal zero at many interest rates. 


1] 
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Figure 5.--Representative types of present net worth functions. 


if the analyst initially sets the interest rate increment large enough, the 
program stop may be avoided and convergence will occur. 


Type 9.--Present net worth is zero at two interest rates. IROR will locate 
only point a. 


Type 10.--Points will be iterated until point b is reached, at which time 
the following message is printed and the program proceeds to the next alternative. 


PNW AS A FUNCTION OF INTEREST RATE IS NOT TENDING TOWARD ZERO 
WITHIN A PORTION OF ITS RANGE FOR ALTERNATIVE X 


The final output is a table listing interest rates and present net worths 
for each alternative until the internal rate of return is reached (see fig. 6). 
When the internal rate of return is found, the precision level applicable to 
that rate is indicated in the present net worth (PNW) column. 


DETERMINING RATE OF RETURN== RESEARCH PAPER SAMPLE PROBLEM 
PRCBLEM NOow 1 TEST DATA 
PRESENT NET WORTH (PNW) AT ALTERNATIVE RATES OF INTEREST 


MANAGEMENT ALTERNATIVE 


ONE TWO THREE FOUR FIVE 

* RATE PNW * * RATE PNW * * RATE PAh * * RATE PNW * * RATE PNW * 
Oo 2e25COE 02 Qs Zo 702SE O07 Oo 4e5500E 02 De 1.1550E 03 Qo 923509E 92 
025C 309320E 92 0,50 2,2C1SE 07 0.5) 302838E 02 0059 1o.1709E 63 0.59 To1O018E v2 
1.00 208318E 02 120C 1,7654E Q7 1,00 203221E 02 Os Oo 1099 523560E 02 
1.250 2,6477E OZ 1.5¢ 1.27e51€ O7 1,50 1.5883E 02 O> Qs 1259 3eS906E 02 
2200 2o47T84E C2 20CC 1.C53SE 07 2000 1e¢0265E 02 Oo Oo 2009 2e9154E U2 
2050 223223€ 92 2050 ToESTIE 0€ 2050 FeS556E O1 Oo Oo 2059 200631E U2 
320¢ 201784E 02 2200 5,14SSE 0€ 3509 226496E O1 o> Ve 3299 123831E 02 
3250 200456E 02 7,50 2eS7T0EE OF 3050 1261855 00 0» Je 3050 @e3747E Ul 
4200 1,9Z28E 92 4000 1.C784E 06 3053 1eCQQ0E 00 Oo Os 4,09 3.9720E O1 
4050 1,8C92E 92 4032 1,CCOCE Cl Co Co or Oo 4059 4e0187E U0 
5200 1.7040E 02 Co Co Os c %o Os 4,56 1,0000E 00 
5050 1,60€5E 92 Ue Qs Co Co Oo 0. Oo De 

6200 1e515SE 02 Co Qe 9 Co Do Oo Co ) 

6050 1.4316E 02 Co o Co Ce Oe Oo Ue Oo 

T2000 1.3536E 02 Cc Oo Co To Os 9 C Qe 

7250 1,2808E 02 Co Co 0. Co Oo 06 Oo 0 

8.00 1,2129E 02 Co c Co Oe De 0 c 

8e50 1,149€E 92 Co O. Co Co i} c Le 8) 

S200 1e09C4E 02 Oo Co 0 c O» 0 Co u 

9250 1.0351E 02 0. Oe Co Ce Oo Oo Ue i) 

10.00 S$e8338E C1 Co Oo C5 Ce O60 Os Co 0) 

ENO OF RUN 


Figure 6.--Fourth report of the IVST computer program showing results 
of use of the tnternal rate-of-return eriterton. 


LIMITATIONS OF IVST 


There are several limitations which must be observed in the use of IVST: 


1. The number of alternatives analyzed within a problem must be less than 
or equal to 20. 


2. The type of series (either terminable or perpetual) must be the same 
for all alternatives within a problem. 


3. The number of years in the investment series for any alternative must 
be equal to or less than 140. 
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4. The length of period to which periodic annual costs or returns apply 
must be identical for all alternatives within a given problem. However, the 
number of periods may vary from one alternative to another. The number of years 
in each alternative (i.e., the investment series length in years) is therefore 
equal to the product of number of periods and period length. 


5. The range and increment in interest rate must be identical for all alter- 
natives within a given problem. 


6. The beginning interest rate must be greater than -100 percent. 


7. The ending interest rate must be less than 85 percent if the investment 
series length for any alternative in the problem is 140 years. Higher interest 
rates may be analyzed if the maximum investment series length decreases. Users 
should consult with computing center personnel regarding the limits on the expo- 
nential function subroutine being used, as this is the factor that limits the 
range of interest rates which may be analyzed. 


8. The interest rate range and increment selected by the user must not 
result in more than 200 interest rates being generated. However, more than 200 
interest rates may be generated during the bisection procedure, which is not 
under the user's control but is automatically controlled by the IROR subroutine. 


9. An attempt to analyze a zero interest rate for a perpetual series will 
result in the program skipping that rate and passing on to the next interest 
rate within the range selected by the user. No program stop or computational 
error will result. 


OPERATION OF IVST 


Input card formats, program flow charts, and a listing of the IVST source 
program are included in the Appendix. IVST is coded in FORTRAN IV and is oper- 
able on either the IBM 7040 or the IBM 360/50 (both time-sharing and nontime- 
sharing modes) computing systems. Slight modification may be necessary before 
IVST can be used with other systems. A total of 31,988 words of core storage 
are required on the IBM 7040 computer for this program. Use of IVST on com- 
puters with smaller memory capacities will require modification of the program. 
The sample problem shown in this paper required |] minute and 4 seconds of execu- 
tion time on the IBM 7040. 
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APPENDIX 


Preparation of IVST Input Cards 


There are nine types of input cards required by the IVST program. The 
format of these input cards is indicated in table 1. Card columns for which no 
punch is indicated should be left blank. All data except alphanumeric must be 
right-justified (i.e., the units position must be at the extreme right of the 
field). It should be pointed out that input card types 1-6 and 9 contain iden- 
tification and control data, whereas input card types 7 (periodic annual costs) 
and 8 (periodic annual returns) are sets of data applicable to independent alter- 
natives. All costs are read for a given alternative. Then all returns are read 
for that same alternative. This sequence (i.e., all costs and then all returns 
for each alternative) is repeated until all alternatives have been read before 
the next input card type is read. The last card contains the end-of-problem 
indicator (MEND). 


BCA Optional Output 


The discounted returns and discounted costs at each interest rate may be 
printed if the user desires. The only modification required in the BCA sub- 
routine to provide for this output is to reproduce BCA program cards BCA 0370 
and BCA 0375 (see listing of the BCA subroutine) deleting the "C" in the first 
column. 
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wu 


Input card 
type 
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] 
2 


STD(I) 
BRATE 
FRATE 
DELTA 
LX 
INIT 


LI(I,L) 


KCXX 


COST(L,KC) 


REV(L,KC) 


MEND 


Table 1.-- \VST Program Input Format 


Alphanumeric study identification. NS seeps A 
Beginning interest rate, decimal. XXX 
Final interest rate, decimal. XXX 
Interest rate increment, decimal. XXX 
Number of alternatives in problem. XX 
Indicator of an initial cost or return; = 0, if 

there is an initial cost or return, = 1, if not. x 
Alphanumeric identification of alternative L. AAAAAAAA 
A second card is needed to handle over 10 alter- 

natives. 

Problem number. XXXX 


Type of calculation; = 01, if series is perpetual, 


= 02, if series is terminable. XX 


Investment criterion; = 01, internal rate of return, 
= 02, benefit/cost and present net worth, = 03, all 
three criteria. XX 


Length of period, years. XX 
Maximu., number of years in the investment series of 
any alternative. If initial costs and returns are 
read, present year should be included in KCXX. XXX 


Number of years in the investment series of each 
alternative L. If initial costs and returns are 
read, present year should be included in KCX(L). XXX 


Alphanumeric problem identification. Avrreteyeytexe A 


Periodic annual cost of alternative L, for period 


KC, dollars. 2 
Aad Dawal CARDS ARE XXXXXXXXXX 


h ; 5 
eG Vited +e ead gare Fok PMOoRE Than % Perieds, 
Periodic annual return of alternative L for period 

KC, dollars. 


Qa PECTS . ened HT 
5 Py / 
handJerlixepatvernattyes: Add)'t/0We/ cakes BYe XXXXXXXXXK 
PeQuires To fea ote Fey More Than fretiods , 
erminal card code; = 98, to do another problem. 

= 99, end-of-run. XX 


Card 


columns 


18-20 


1-10, 11-20, .. 


1-10, 11-20, .. 


2-3 


.-, 598-60 


-» 71-80 


-, 71-80 


FLOW CHART FOR MAIN ROUTINE IVST 


LEGEND FOR FLOW CHARTS 


START OF DO LOOP i 
TOSS! (SRR 


TERMINUS OF LOOP i ess PROGRAM 
BEGINNING AND END 


INPUT OR OUTPUT DIRECTION OF FLOW 


DECISION 


MPUTATION AN ENTRY FROM, OR EXIT 
oe PROCESSING TO, ANOTHER PART OF 
THE PROGRAM FLOW CHART 


Page 1 of 7 


DIMENSION VAR/ABLES 
INPUT FORMATS 
ovTPuT FORMATS 


KCT=KCxx / IP 


2] 4 


Page 2 of 7 


COST(L,KC)=0 
REV(U,k¢e)=0 


BRATE FOO 
AND 
FRATE #0.0 


BRATE= BRATE+ DELTA 
BRATE 
VES 
FRATE= FRATE— DELTA 
YES 


NR= ((FRATE-BRATE)/ DELTA) |-08 
RTLOG(1) = ALOG (1.0 + RATE (i)) 


4 


Page 3 of 7 


ACSA) NR 


RATE (T)= RATE (T-1)+DELTA 
RATE (x)= RATE(T) + DELTA a 


RTLOG (Z)= ALDG (1.0+ RATE (1S) 


READ ACI) 
as =1,/3 


Sly 


KeS=KC\XW)/IP 


READ CoST(L, KC) 
REv (L,k¢) 
KC=1,KCS 


Pe: 


23 


24 


eA) Page 4 of 7 


WRITE 
PERPETUAL 


Sale 


Page 5 of 7 


LI=LI+5, | baka Wome lee ee a 
L3=L3+S, L4=14 +5 — 


L5s=L545 


No By 


NES 


eA) 25 


Z 
: 
2 
C 


26 


BECO GEiEe)— Cosme Gane 
RETURNCL Le)= REV (4, kc) 
YR (14,42) = LE 


LXX= Lx X+] 


> > 
> 


Page 6 of 7 


YRCLL2) = YR (L, Le)- 1.0 


Page Or i 
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FLOW CHART FOR SUBROUTINE BCA 


Page 1 of 5 


DIMENSION VARIABLES 
COMMON STATEMENT 
FoRMATS 


a LLL =I, LX 


ILTI=A\,nR 
VALIN (LLL,IIT) =0.0 


ie BCR (421, 7z25r=0,0 
<a ele 


LX 


Kex A= K0Y CL) 
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De =DKC+PECO (L,K0)/EXP(RTLOG(L) YR (i, KC)) 
DKI = BKI+RETURN(L, ke) /EKP(RTLo GZ)¥ YR (1, ke)) 


DKC = DKC+ PEC (L, ke) 
DK|= DK/ + KETU RN (Ly ke) 


DISCO= EXP(RTLOG (I)* KLY) 
DIS!= Gox |.O/ DISCO =1.0)) 
Oki= DK) xDIS 
Dke= DKC¥DIS 


VALIN (4,1)= DKI - OxKe 
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By Page 3 of 5 


BCR(L, I) =DKI/DKG 


I=|,NVR 


3] * 


32 


LINE=LINE+] 


cy 


yes 
No 


Page 4 of 5 


Page 5 of 5 


LizLi4S° ,La=La+sS) 
La=L35 ) L4=J45 5 
LS=L54+5 
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FLOW CHART FOR SUBROUTINE IROR 
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Sw 


DIMEN SION) VARIABLES 
COMMON STATEMENT 
FoRMaTS 


(3 \ 


ae 
LULe I) LX 


ITI=nNerR 


KCXA= KCX (L) 
XLY= FLOAT (Kcx (L)) 
= 2 
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T=!, NR 


RTLOG CZ) = 0.0 58) 
YES 


NO 
K@=I, KC XA 


OKC=DKC+ PECOCL,KC) /EX PCRTLOG (zk YR Lk) 
| OK I= DKI+ RETURNGL, ke) /EXP(RTIDG (Z)* YR(L, ke)) 


> 


TBVAL= DbkK]—DKC 


VAUN(L,I)=TDVAL 
YES 
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DISCO = EXPCRTLOG KXLY 
VALIN(L,Z)=T DVALK(.0+1. 0/(01S¢0 -/.0>) 


VALINL,I) £0.01 
AND 
VALIN(L,I)= ~9. 01 


y 


YES 


\\ 


= Oo 


yy BY YY & 


37 


38 
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DKC =DKC + PECO (L, Kc) 
DK| = DKI + RETURN €L, KC) 


NRATE(L) =) RT 
X RATE (IRT)=FRATE 
LR=| 


LINE= LINE+] 


| 
Vv 
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MRATE(L)=IRT 
XRATE (L,IRT) = RATE (LRT) 


VALIN(Z,TRT) £0. 0| 


AND 
ALIN(L,rRT) 2—24L 


VALIN (L,IRT) =O.0 


4] % 


YCI)= vaLin (L,IRT) 
Y¥(4)= VALINCL ,ERT-1) 


XC1)= RATE (LRT) 
: IT=3, /00 


X (a) = RATE (ZRT-1) 


Yes 
NO 


X@z)=( @I-a) +x(r1-1)) 2.0 
RLOG= ALOG (1.0+X (TB) 
DISCO= EXP (RLOG * XLY) 
DKC = 0.0 
DKI = 0.0 


— 


Kc=| p 


OKC= DKC+PECOC,ke)/ EXP (RLOG¥VR(L,kC)) 
DK i= DKI4+RETURN(L, KC) /EX P(RLOG-& YR(L, K0)) 


wb 
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TOVAL= DKI- DKc 


- 


Y(II)=TAVALK().0+/,0 /bIsco-|.0 


YI) < Py 
AND 
y (TI) > PL VES 


x (Ir-!)=x GT-a) IS 


43 
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& 
€2,—_D 
= 


VALIN(L, IRT) = (ID) 
PH= PH 10.0 
Pi= PES Je.0 


YES 


ws 


A Page. 11 of 14 
¥ RATE (L,IRT)= xD) 


{3 


NRT = NRATE(L)— | 


IM= 1, NRT 


KK=NRATE(L) — TIM 


oS YRATE (Ly KK)= RATE(KK) 


€ 
ee 8 


No 
GA > 
YES 


KM=LX -! GS 


oS EKaNs | 
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SHE lo Say 


ITEM P=MRATE (71) 
NRATE (T1)=NRATE(T Hi) 
NRATE@Ie) = LTEMP 


NRR=NRATEC 


(tou 


Ws 


=S 
=4 
aS 


I=I, VRR 
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IVST SOURCE PROGRAM 


GLIOLSAI$7 SATLVNYUIL IV NI LINGOYd JO GISIA ANIWA IANNY DJIGOTYAd=(9'1) A4x 
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OLTOISAI 
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OSTOISAI 
SYTOLSAT 
OYTOLSAT 
SETOLSAI 
O€TOISAT 
SZTOLSAI 
OCTOLSAI 
STTOISAI 
OTTOLSAT 
SOTOLSAI 
OOTOILSATI 
G6OOLSAI 
O600LSAI 
68001SAI 
B8O00LSAT 
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Headquarters for the PACIFIC NORTHWEST FOREST AND 
RANGE EXPERIMENT STATION is in Portland, Oregon. 
The Station’s mission is to provide the scientific knowledge, 
technology, and alternatives for management, use, and 
protection of forest, range, and related environments for 
present and future generations. The area of research encom- 
passes Alaska, Washington, and Oregon, with some projects 
including California, Hawaii, the Western States, or the 
Nation. Project headquarters are at: 


College, Alaska Portland, Oregon 
Juneau, Alaska Roseburg, Oregon 
Bend, Oregon Olympia, Washington 
Corvallis, Oregon Seattle, Washington 


La Grande, Oregon Wenatchee, Washington 


4407) 


The FOREST SERVICE.of the U:S. Department of Agriculture 
is dedicated to the principle.of multiple use management of the 
Nation’s forest resources for*sustained yields of wood, water, 
forage, wildlife, and recreation. Through forestry research, 
cooperation with the States. and. private forest owners, and 
management of the National. Forests and National Grasslands, it 
strives — as directed by..Congress —<to provide increasingly 
greater service to a growing Nation. 
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